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The quantum 0{2) model in 2 + 1 dimensions is studied by simulating the 3d 0(2) model near 
criticality. Finite densities are introduced by a non-zero chemical potential fi, and the worm algo- 
rithm is used to circumvent the sign problem. The renormalisation is discussed in some detail. We 
find that the onset value of the chemical potential coincides with the mass gap. The /^-dependence 
of the density rules out Bose-Einstein condensation and might be compatible with an interacting 
Fermi gas. The fi-T phase diagram is explored using the density and the magnetic susceptibility. 
In the cold, but dense regime of the phase diagram, we find a superfiuid phase. 
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I. INTRODUCTION 

The 0(2)-model as a statistical field theory in three di- 
mensions has always enjoyed great phenomenological im- 
portance: to name only a few applications, it illuminates 
the superfiuid transition (A -point) of pure '^He [1], it de- 
scribes the electromagnetic penetration depth in certain 
high- Tc superconducting materials [5] and it might pro- 
vide insights in the Bose-Einstein condensation of atomic 
vapour (hard sphere bosons) [3]. It is well known (see 
e.g. [1]) that the 3d classical theory possesses a sec- 
ond order phase transition at a critical inverse temper- 
ature /3c = 0.45420(2) separating the disordered phase 
at /? < /3c from an ordered state for /3 > /3c. In solid 
state physics, the atomic characteristic length a sets the 
fundamental scale of the theory. For later use, we intro- 
duce the dimensionless correlation length as the physi- 
cal correlation length in units of the atomic scale, i.e., 
Ciatt(/3) ■— £,/a. It is this correlation length which di- 
verges at the phase transition. 

The 0(2)-model inherits its name from a global 0{2) 
symmetry, which lets us define a (Noether) current. The 
generalisation of the 0(2)-model to include a finite chem- 
ical potential, which favours the proliferation of density, 
has recently attracted a lot of interest: The probabilistic 
weight is now complex, and the model inhibits the nu- 
merical study with classical Monte-Carlo methods based 
upon importance sampling with respect to a real and pos- 
itive probabilistic weight. Thus, the model shares this 
notorious sign problem with the theory of strong inter- 
actions, i.e., QCD, at finite baryon densities. The QCD 
phase diagram as a function of temperature and chem- 
ical potential is still not known form first principle cal- 
culations, and the prospects are bleak unless a change of 
paradigm let us address theories with sign problems on 
more generic grounds. Such breakthroughs might equally 
well be made in the context of simplistic theories such 
as the 0(2)-model at finite densities. Over the recent 
past, quite some progress has been made in this direc- 
tion: Monte-Carlo simulations which are based upon im- 
portance sampling with respect to the density of states, 
an always real and positive quantity, might circumvent 



the sign problem O [6] . The so-called complex Langevin 
approach [71 18] is based upon a complexification of the 
fields and might be ideally suited to address complex ac- 
tion systems. This approach has been largely explored 
over the recent years [9j [10] . Concerns about its relia- 
bility have been raised recently, since the approach, al- 
though it converges, it does not give the correct answer 
in certain cases. While the Bose gas is a spectacular suc- 
cess story for the complex Langevin approach fill, 112], the 
approach fails to produce the correct answer for a close 
relative, the 0(2)-model in certain regions of parameter 
space |13j . Progress has been made recently by devis- 
ing criteria for correctness of the method 14 . Another 
promising idea, firstly put forward by Chadrasekharan 
in [TS], is based upon a reformulation of the theory in 
terms of new, potentially non-local variables in a sign- 
problem free manner. In particular, the so-called world- 
line or worm algorithm ,16, il7j has been identified as an 
approach with a wide range of applications including the- 
ories with sign problems. As pointed out by Gattringer, 
the so-called SU{3) spin model, which is motivated from 
dense QCD in the limit of heavy quarks, is accessible by 
means of worm (also called flux type) algorithms even 
at finite densities [T8Vl21j . Most relevant for the study 
here is the flux representation of the 0(2)-model at finite 
densities. A thorough analysis of the classical model has 
been presented by Chadrasekharan in [22, . 

Although a range of classical theories with sign- 
problems have been explored in the recent past, to the 
best of our knowledge such an investigation of a finite 
density quantum field theory has not yet been performed. 
In the present paper, we will bridge this gap and system- 
atically explore the quantum limit of the 0(2)-model at 
finite densities. For this purpose, we rely to a large extent 
on the fiux type algorithm [52] which will give us access 
to the full phase diagram as a function of temperature 
and chemical potential. Generic quantum phenomena 
such as superfiuidity will be described by a first principle 
simulation. 

On a practical note, the quantum limit of the model 
is accessible using simulations of the classical theory and 
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by performing a detailed scaling analysis near the critical 
point. This will then reveal the properties of the quan- 
tum 0(2) theory in 2 + 1 dimensions. Let us explore the 
correspondence between classical and quantum theory a 
bit further: In the latter context, the inverse temperature 
P is re-interpreted as the spin-spin coupling strength, and 
the temperature T is now attained by the extent of the 
torus in "time" direction. The inverse correlation length 
is now interpreted as the mass gap, to = 1/^. This mass 
takes over the role of the fundamental scale of the quan- 
tum field theory and is kept fixed under a change of the 
coupling strength /3. The formerly introduced atomic 
length a (also called the lattice spacing in the quantum 
context) merely plays the role of a regulator. This regu- 
lator a now necessarily depends on the coupling strength: 
ma{f3) = l/^iatt(/5)- Tuning the coupling /3 to the criti- 
cal value (at a fixed value for to) implies that the lattice 
spacing tends to zero. This then installs the continuum 
limit of the quantum theory. 



II. MODEL SET-UP AND SYMMETRIES 

In 2-1-1 dimensions, the 0(2)-model, also called the xy- 
model, is formulated on a cubic grid of size V = NtN^Ny. 
Dynamical degrees of freedom are the angles (j)x which 
are associated with the sites x of the lattice and which 
parametrise the unit vectors 



Action and partition function are given by 
So= ^ nlnx+„ = ^ cos(( 

t=(xu) t=(xu) 



(1) 



4>x+iy) , (2) 



Z = jvct>, expl/JS-ol^]}, 



(3) 



where /3 is the coupling constant, and where x -I- is 
the neighbouring site of a; in direction. The sum in 
([2]) extends over all links £ of the lattice. The model 
possesses a global 0{2) symmetry. 



(4) 



where 0((/>) G 0(2). This displacement symmetry gives 
rise to the Noether current: 



Indeed, defining the divergence of the current by 
A^j^(a;) oc ^{x) - j ^{x ~ fi) , 



(5) 



(6) 



and using the identity (for any xq) 



'D<l>x 



d 



exp 



0, 



we find that any number of insertions of the divergence 
^ into the partition function causes the integral to van- 
ish, e.g.. 



A^j^(a;o) exp|/3 5o[0]| 







(7) 



The model can then be readily generalised to finite charge 
densities, by introducing a chemical potential /i to the 
action in the usual way |23j : 

S[<t>] = X! ^^^i't^x - ipx+v - Ut-SpQ) ■ (8) 

t=(x,v) 

Thus, the imaginary part of the action is proportional to 
the conserved charge ja{x) ([5]), i.e.. 



S[cf>\ = 



E 

= {x,u] 



|^COsh(/x) COii{(t>x - (t>x+y) 



(9) 



i sinh(/i) (5^0 sin( 



(jix+v) 



For large values of the coupling /3, the global 0{2) sym- 
metry Q is spontaneously broken in the infinite vol- 
ume limit. The operator coa{(j)x) is not invariant under 
the symmetry transformation, and its expectation value 
serves as an order parameter in the infinite volume limit. 
Note, however, that for any finite volume V , the expec- 
tation value 



0, 



(10) 



i.e., it vanishes for any value of the coupling /3. This can 
be easily seen by performing the variable transformation 

cj)x^(t>x+V (11) 

in the integral (|3|. In order to trace out any indications 
for spontaneous symmetry breaking at finite volumes, we 
add a source term to the action which breaks the 0(2) 
symmetry explicitly. We then study the response of the 
expectation value in ( 10 1 to variations of the external 



source. Here, we choose 

Z[j] = I Vc^x exp{/?5[0] + J ^cos(0,) }. (12) 

X 

The response function is defined by 

R,{P) = (cos(0,)) = In^b'] . (13) 

To decide whether the 0(2) symmetry is spontaneously 
broken, we take the limits and find 

,. ,. r. f = Wigner-Weyl realisation 
hm hm i?, (/3) <^ , „ ^ , u 1 • 

v_>oo I ^ spontaneous sym. breakmg , 

where the order of the limits is crucial. For this program, 
it is sufficient to study small perturbations j only, which 
leads us in leading order and finite volumes to 

R{3) ^ xiP) J + 0{f) , (14) 
= ^^lnZ[j]|j=o = {^cos((/):ro)E'^°^('^^)/' 
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FIG. 1. The average action as a function of /3 for the 3d 
0(2) model using a 32'^ lattice. 



FIG. 2. The magnetic susceptibility per volume as a function 
of /3 for the 3d 0(2) model using a 32'^ lattice. 



where x is (half of) the so-called magnetic susceptibility. 
If X remains finite in the infinite volume limit, the re- 
sponse function vanishes in the limit of vanishing source 
j, and the system is in the Wigner-Weyl phase. Hence, 
a diverging magnetic susceptibility in the infinite volume 
limit is a necessary condition for the spontaneous sym- 
metry breakdown. 

Note also that x can be viewed as the expectation value 
for J2 cos(0) where the factor cos{(l)xo) acts as a reference 
against which the spin orientation is counted. The mag- 
netic susceptibility can be written in a manifestly 0(2) 
invariant form. To this end, we rewrite the expectation 



value on the right hand side of ( 14 ) as 



in 



The latter identity explains the role of x(/3) as an or- 
der parameter for the spontaneous breaking of the 0(2) 
symmetry: in the disordered phase, the spin correlation 
function exponentially decreases over the distance of the 
correlation length ^, and x{P) is independent of the sys- 
tem size. If the correlation length near criticality exceeds 
the system size, we find by means of the sum at the right 
hand side of (17) that x[P) diverges with the volume. 



III. THE STATISTICAL FIELD THEORY 



A. Algorithms and first results 



Wolff cluster algorithm 



Using the variable transform (11), we conclude that the 



latter expectation value vanishes for j — Q (at any finite 
volume). Thus, we can write the magnetic susceptibility 
as 

X 

which is manifest 0(2) invariant. Introducing the spin 
correlation function, by 



G{x ~y) = (rix-ny^ , 



(16) 



the magnetic susceptibility can be easily related to the 
integrated spin correlation function: 



(17) 



For vanishing chemical potential, i.e., /i = 0, and van- 
ishing source j, the Wolff cluster algorithm [24] offers a 
very efficient way to simulate the 0(2)-model. The rea- 
son is that the clusters, which are updated in one Wolff 
step, grasp the physics of the model. This implies that 
the auto-correlation time depends only weakly on the sys- 
tem size as refiected by a dynamical critical exponent 
close to zero. The algorithms proceeds as follows: 

1. Choose a random unit vector r; 

2. Activate the link between two neighbouring sites x 
and y with probability 

P{x, y) = I - exp {min [0, -2 l3 {ri'^ ■ f) {uy • r) ]} ; 

3. Select a random site Xq and find the subset C of all 
spins which are connected to (pxa by activated links 
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FIG. 3. Detail of the magnetic susceptibility as a function 
of /3 close but smaller than the critical value 



4. Replace 

— > ria; — 2 (tT^ ■ r)r V a; G C . 

Observables are best calculated using the so-called im- 
proved estimators |25l 126 ) . In particular, the magnetic 
susceptibility can be directly related to the cluster prop- 
erties and as such calculated "on the fly" : 



(18) 



x(^C 



where C denotes a particular WolfF-cluster, and Vc is the 
number of sites belonging to this cluster C . 



2. 1-shot heat bath algorithm 

Heat bath algorithms offer an easy access to expecta- 
tion values formulated in the terms of the original spin 
degrees of freedom. We here use the heat bath algorithm 
to test and benchmark the flux algorithm outlined in the 
next subsection. For this purpose, we have developed a 
heat bath algorithm with 100% acceptance rate, which 
has been inspired by the method from Bazavov and Berg 
in P?'. In step 1, we choose the spin (t)^ for the local 
update step. If {xy) denotes the links joining the sites x 
and y, the action can be written as 



\P ^ cos{4)y) + j cos{(j)^) + 

ye{xy) 



[f3 Yl sin(0,)' 



y<£{xy) 



A cos{(l)x - -0), 



where A and ip depend on /3, j and the neighbouring 
spins. In step 2, we define 



exp 



I A cos{v — "0)1 dv 



We then choose a random number u E [0, 1] and solve 
the equation 

for the new value (j^x for the spin at site x. In practice, 
the integration was performed using the Gauss Legendre 
method, while the non-linear equation was solved using 
bisection. 



3. Worm algorithm 

For finite values of the chemical potential, i.e., /i 7^ 0, 
the action ([s]) possesses imaginary parts. The theory is 
hampered by the infamous sign problem and cannot be 
simulated at reasonably sized lattices using the cluster or 
the heat bath algorithm. The 0(2)-model is one of the 
rare cases where simulations are nevertheless feasible by 
using the so-called worm algorithm [H [T71 [H [55] • The 
algorithm has been thoroughly discussed in [52] . We here 
briefly review the algorithm for the 0(2)-model at finite 
chemical potential before we take it a step further to 
derive the Kramers- Wannier duality in the next section. 
Using the identity 



cos(0) 



E 

fc— — 00 



hW) e^'* , 



(19) 



where h- = I-i- is the modified Bessel function of the 



first kind, the partition function ( 12 1 becomes 



Z[0]=y" Vcjjx e'^'^1'^1 =7^*^ 

{ke}l={xy) 

{ki}e={xy) \iex ) 



(20) 



where the notation £ e x refers to all links which are 
attached to the site x. If an integer "flux" ki is associated 
with each link I of the lattice, the (5-function constraint 



in ( 20 ) implies that the total (directed) flux which enters 



a site must vanish. 

Including the case of a finite chemical potential /x, the 
worm algorithm has been widely used in [22] to explore 
properties of the 0(2)-model. The algorithm works as 
follows 22]: 

1. Pick a starting point x and set the counter c = 
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FIG. 4. The percolation probability Pporcoi ( 25 1 as a function 
of /3. 



2. Choose at random one of the 6 directions (includ- 
ing negative ones); let y — x + v he the neigh- 
bouring site associated with the chosen direction 
{h' = ±1, ±2, ±3) and £ the associate link (xy) = £. 

3. If a positive direction was selected, with probability 
exp{nS^o)Ikf,+i{l3) / IkdP) change kg to ki + I and 
move to the neighbouring site y. Otherwise stay. 
In any case, increase c ^ c + 1. 

If a negative direction was selected, with probabil- 
ity exp{-iJ,S„o)Ike-i{l3)/Iki{l3) change ke to ke - 1 
and move to the neighbouring site y. Otherwise 
stay. In any case, increase c ^ c + 1. 

4. Stop if the current position coincides with starting 
position, i.e., it y — Xg- 



The average action is easily obtained by 



d_ 

9/3 



E 



(21) 



It has been argued in [2^, that the magnetic suscepti- 
bility is given in terms of the average number of counts 



XW) = ^ (c) . (22) 



4- Action average and magnetic susceptibility 

High precision simulations of the model have been 
reported in the literature for zero chemical potential 
e.g. in [551 US and for finite densities in We will 



benchmark the outcome of our algorithms against these 
findings before we will move on to study the quantum 
field theory (QFT) limit at finite densities in the next 
section. We will employ 

• the Wolff algorithm to obtain the correlation length 
in units of the lattice spacing, which will set the 
QFT scale, 

• the worm algorithm for the simulations at finite 
densities, 

• the heat bath algorithm to crosscheck the numerical 
results. 

To begin with, we have compared all three algorithms by 
calculating the average action as a function of the cou- 
pling /3 using a 32'^ lattice. Our findings are shown in 
figure [T] We find a nice agreement for all three algo- 
rithms for a range of /3 values covering the disordered 
phase for /3 < /3c, /3c = 0.45420(2) g^, as well as the 
ordered phase /3 > /3c. Our results for the magnetic sus- 
ceptibility X a-re shown in figures [2] and [3] The transi- 
tion from the Wigner-Weyl phase to the phase with a 
spontaneously broken 0(2) symmetry is clearly visible in 
figure [2j As expected the Wolff and the worm algorithm 
work flawlessly in the broken phase while the results from 
the heat bath algorithm are plagued by large autocorre- 
lation times. The region of /3 just a bit smaller than /3c is 
of particular interest for the QFT limit. Here, all three 
algorithms perform well (see figure [s]). A good agreement 
with the fit 



Xo(l-£) , P^Pc, (23) 

Xo = 0.5045(1), /3c = 0.45420(2), 7 = 1.324(1), 



x(/3) 



presented in [U [25] is observed (note, there is a factor of 
1 /2 difference in the definition of our magnetic suscepti- 
bility and that in [ill^). 



B. Winding sectors and dual theory 

1. Flux percolation 

By virtue of the periodic boundary conditions, the 
closed flux lines of the worm algorithm might wind 
around the torus. Since we currently work with a sym- 
metric lattice of equal size N in all directions, we con- 
fine ourselves to the study the winding in time direction. 
The winding number n„can be calculated during the con- 
struction of particular fiux line. Here, we only need to 
keep track how often the fiux line leaves the torus at the 
Nt = N and Nt = 1, respectively. Every time the flux 
line leaves the torus d.t Nt = N in positive time direction, 
we increase riw by one unit while if the flux line leaves the 
torus at Nt = 1 negative time direction, we decrease riw 
by one unit. One naively might expect that the average 
winding number is conclusive on the winding of the flux 
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FIG. 5. action density for the standard 3d-0(2) model in 
comparison with model with winding number n^n ~ only. 



lines. This is, however, not true due to the abundance of 
small size flux lines which all carry riw = 0. The more 
interesting quantity is the percolation probability Pporcoi: 
if we choose randomly a particular link of a flux line con- 
figuration, how big is the probability that this link is a 
part of a flux line which at least winds once through the 
torus? If L denotes the number of links of a particular 
flux line and if 



L for 
else 



^0 



the percolation probability is given by 
-^pcicol — {Lu.) I {L) . 



(24) 



(25) 



Figure [4] shows our numerical result as a function of j3. 
Note that due to the integration of the 0^; in ( 20 1 all ref- 



erence to the global 0{2) symmetry has been lost in the 
flux line representation of the partition function. Here, 
We find empirically that the flux lines do not percolate in 
the Wigner-Weyl phase while the flux lines starting wind- 
ing around the torus in the phase with spontaneously bro- 
ken symmetry. This is intuitively clear by noting that the 
correlation length diverges for /3 — ;> /3c. In the flux line 
representation, this long range correlation can only be 
provided by the "infinite" extent of the flux line clusters. 



2. Kramers- Wannier duality 

The reformulation of classical field theory in terms of 
flux variables has been dubbed the dual formulation in 



the literature (see e.g. [30]). By its nature, this refor- 
mulation is non-local, using the original lattice, and it 
is therefore quite different from the generic Kramers- 
Wannier duality, which employs local degrees of freedom 
related to the dual lattice. Starting from the worm al- 
gorithm as formulated in subsection [III A"3l we here de- 
rive the dual formulation of the dense 0{2) model in the 
Kramers- Wannier sense. We will find a close relation to 
the flux formulation: the dual theory turns out to be 
a gauge theory with the fluxes being the non-local but 
gauge invariant degrees of freedom of this formulation. 



Oi 



Starting point is the flux formulation (20). 
will be to resolve the (5-function constraint on the flux 
elements ki using local variables. Naturally, the emerging 
dual theory is free of the sign problem at finite chemical 
potential as the flux line formulation (20) is. To this end, 
each element of the lattice, i.e., site x, link plaquette 
p and elementary cube c, is mapped in the usual way to 
the dual lattice: 



P, 



c — ^ X 



where the "bold" quantities denote the entities on the 
dual lattice. Fields depending on the entities are mapped 
as well, e.g., 

fcf fcn . 



The 5-function constraint in (20) then becomes on the 
dual lattice: 



^h = ^ ^fcp-O, 



(26) 



where p S c addresses all plaquettes p which build up 
the faces of the cube c. The crucial point is that this 
constraint can be solved locally for the theory with zero 
windings, fcp are plaquette variables on the dual lat- 



tice, and the latter constraint (26) is reminiscent of the 



Bianchi identity of a Z-gauge theory. Hence, we solve the 
constraint by introducing the integer valued link field Zg 
of the dual lattice by 



(27) 



The definition of the links which represent the pla- 
quette fcp is not unique. The gauge transformed links. 



Z 



where ri(x) is integer valued, still give rise to the same 
plaquette in (27). A Z-gauge symmetry arises in the link 



representation of the dual theory. The partition function 
(20) of the dual theory can be written as: 



{Zi} P 



fc„ expi ^^^(pgspatial) ^p} • (28) 



The dual formulation sheds light onto the "worms" of 
the flux line representation. By virtue of (27), this flux 
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is given by the plaquettes of the dual Z-gauge theory. On 
the other hand, the dual plaquette defines the vorticity 
of the dual theory. Vorticity is conserved by the Bianchi 
identity, and non-trivial vorticity conies form closed sta- 
ples of plaquettes. Hence, we find that the closed flux 
lines, i.e., the "worms", appear to be the vortices of the 
dual formulation. 

So far, we have established a local equivalence between 
the standard 0{2) model and the Z-gauge theory on the 
dual lattice. For a complete correspondence, we also have 
to match the boundary conditions of both formulations. 
This a non-trivial task: let us consider, e.g., periodic 
boundary conditions for the dual theory. If A is the max- 
imal planar surface in the xy-plane and (A) its boundary, 
we find that the total flux through this plane necessarily 
vanishes: 

Yl Ze = 0. 

On the other hand, the standard formulation of the 0{2)- 
model with periodic boundary conditions does not re- 
strict the worm configurations to the sector of zero wind- 
ing in z-direction. Obviously, we have to give up periodic 
boundary conditions for the Z-gauge theory if we want 
correspondence. If is the worm winding number, the 
corresponding Z-gauge theory simulating this sector fea- 
tures non-trivial boundary conditions such as 

Z2{x = 1, 1, 1) = Z2{x = N,, 1, 1) + n^, 
(periodic else). 

The issue, however, is that it is a priori unknown with 
which weight these winding sectors must contribute in 
order to correspond to the standard 0{2) model with 
periodic boundary conditions. We expect, however, that 
the non-trivial winding sectors are not important in the 
disordered phase. In this case, we would expect a corre- 
spondence between the standard and the dual formula- 
tion both with periodic boundary conditions. 

We have verified this by a direct simulation. The action 
density calculated with the Z-gauge theory with periodic 
boundary conditions (i.e., zero winding number = 
0) is contrasted to the result from the 0{2) model with 
periodic boundary conditions in figure [5j As expected, 
both action densities agree in the disordered phase and 
are different in the broken phase when windings have a 
role to play. 

We find that the 3d 0(2) (or xy) model is dual to a 
theory of vortices. This might explain part of its phe- 
nomenological success for solid state systems the ther- 
modynamics of which are influenced by vortex dynamics 
(see e.g. [1]). 
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FIG. 6. Magnetic susceptibility x in the QFT limit of the 
3d 0(2) model as a function of the temperature T in units of 
the mass gap. 



IV. QUANTUM FIELD THEORY 
WIGNER-WEYL PHASE 



A. Finite temperatures 

Here, we study the emerging quantum field theory in 
the critical limit of the 0(2) statistical field theory. We 
will consider finite temperatures, but we will confine our- 
selves to zero density. 

If we approach the critical coupling /3c from below, 
i.e., /3 S Pc, the 0(2) theory remains in the Wigner- 
Weyl phase, while the correlation length ^(/3) in units of 
the lattice spacing a diverges. For /3 close to the critical 
coupling, one finds a scaling relation of the type: 



m = ^0 



1 



I 



(29) 



The correlation length is best calculated using an im- 
proved estimator and the Wolff algorithm. In fact its 
was shown that the spin correlation function (161 can be 
estimated bv PSII^ : 

G{x^y) - 2/ 



{r-n^){r-riy)e{C,x)6{C,y)) , (30) 



where 



B{C,x) = 



1 for a; e O 
else. 



The scaling parameters have been obtained to good accu- 
racy (see e.g. [4j) by fitting the scaling relation to Monte- 
Carlo data: 



/3c = 0.454157(14), 
^0 = 0.4866(26) . 



0.6711(16), 



(31) 
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FIG. 7. The charge density p as function of the chemical 
potential in units of the mass gap m in the QFT limit. 



For large correlation length, the theory becomes inde- 
pendent from the underlying lattice structure and can be 
effectively described by a quantum field theory (QFT). 
The QFT limit is attained by defining a fundamental 
length scale which plays the role of the only free param- 
eter of this theory. Here, we choose the so-called mass 
gap m of the theory which is provided by the inverse cor- 
relation length. The scaling relation ( 29 1 then implies a 



dependence of the lattice spacing a on the bare coupling: 
ma{(3)^r\(3) ■ (32) 

The "bare coupling /? is tuned towards the critical cou- 
pling. This implements the QFT limit by sending the 
lattice spacing (in physical units) to zero. The bare cou- 
pling f3 is no longer a free parameter of the theory, but its 
role has been taken over by the dimensionful parameter 
m (dimensional transmutation). 

QFT observables are constructed from the "bare" ob- 
servables by scaling the quantity of interest in units of 
the mass gap m (to the power of its canonical dimension) 
and by taking into account the wave function renormali- 
sation. Because of its relation to the correlation function, 
the canonical dimension of the magnetic susceptibility is 
two, and we write: 

ZxiP)xm = Xqfta' , (33) 
where Zy is the wave function (or better composite op- 



erator) renormalisation constant. With the help of (23 1, 
(p9| and (IMl), we find: 



xW) 

em 



Xo 
so 



1 - 



(34) 



In a "minimal subtraction scheme", we would choose 



/3c 



(35) 



where the anomalous dimension for this operator 
turns out to be quite small: 



7 - 2i/ 



-0.018(2) 



(36) 



Needless to say that the actual value of the renormalised 
quantity Xqft must be fixed by a renormalisation condi- 
tion. 

Finite temperatures, i.e., T 7^ 0, are implemented by 
reducing the extent of the lattice in the time direction: 



1 



Nt a(/3) 



(37) 



where Nt is the number of lattice points in time direction. 
Once the observable has been renormalised achieving vac- 
uum values which enjoy a physical interpretation in the 
continuum limit a — >■ 0, we can start to make predictions 
at e.g. finite temperatures. Here, we study the temper- 
ature dependence of the magnetic susceptibility. To this 
aim, we primarily carried out calculations at /? = 0.44 
using a 64^ x Nt lattice. For a study of the lattice spac- 
ing (in-)dependence, we also did simulations utilising a 
32^ X Nt lattice and /3 = 0.414. Note that both simula- 
tions implement roughly the same spatial volume since 

a(/3 = 0.414) « 2a(/3 = 0.44) . 

Our findings are summarised in figure [6] The overall 
agreement between the two curves is satisfactory bear- 
ing in mind that the last points at large T correspond to 
Nt — 2 for which we might expect quite some rotational 
symmetry breaking effects. We observe the magnetic sus- 
ceptibility stays roughly constant for T < 0.4 m while for 
T sa m it has fallen approximately to half its vacuum 
value. 



B. Cold, but dense matter 

Let us now study finite densities by turning on the 
chemical potential. The physical chemical potential /iphys 
is related to bare chemical potential /i, the parameter 
e.g. featuring in (201, by 



A^phys/^ 



(38) 



We will do so at zero temperatures, which are adopted by 
performing simulations using a symmetric lattice, N = 
Nt. Again, we study the QFT limit of the 0(2) theory, 
but now for non-zero values of the chemical potential. 

The first observables which we will study is the matter 
density. The matter density is easily accessible using the 
worm algorithm (see (20l): 



P 



1 d\nZ 
V2 dfi 



1 

V2 



(ko) , 



(39) 



9 



300 



250 



200 



3d xy model, 32 



o 



150 



100 



50 



r 

A 



o 


P=0-45 


□ — □ 


13=0.44 


O 


13=0.43 


A 


13=0.42 




FIG. 8. Magnetic susceptibility x in the QFT limit of the 3d 
0(2) model as a function of the chemical potential /j, in units 
of the mass gap. 
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FIG. 9. Magnetic susceptibility x P^^ total volume as a 
function of the chemical potential /i. 



where V2 = iV^a^ is the "spatial" volume. In physical 
units, we find 



em 

iV2 



(40) 



Since the charge is conserved, there is no wave function 
renormalisation of the density. This is confirmed by a 
direct calculation (see figure [?]) using several (3 values 



ranging from 0.42 to 0.45 implying that the correspond- 
ing lattice spacing a changes by a factor of approximately 
4.2. Yet, the density in physical units nicely falls on 
top of a uniform curve. Most importantly, we observe 
that the density vanishes for chemical potentials smaller 
than the mass gap, i.e., /i < m. With other words, 
the onset chemical potential agrees with the mass gap, 
and the simulation is free of a Silver-Blaze problem |31j . 
Secondly, the density smoothly starts to rise with the 
chemical potential at the onset value. This is markedly 
different from the case of a free Bose gas for which we 
would observe that the density diverges logarithmically 
when approached from above. This would indicate the 
well known phenomenon of Bose-Einstein condensation. 
However, our data are more compatible with a quadrati- 
cally rise of the density with the chemical potential which 
is reminiscent of a free 2d Fermi gas at low temperatures: 



PFo 



Although the standard 0{2) model is formulated in terms 
of bosonic degrees of freedom , it is not ruled out that 
model near the onset transition is well described by an 
effective fermion theory. To trace this out further, we 
have fitted the numerical data in figure [7] to a simple 
scaling law: 



p/m^ = ao(/i/m-l)°i 



(41) 



We find, however, that the data are not well represented 
by this ansatz. Performing independent fits to each set 
for a given /3, the results are collated in the first two rows 
of table H 



fflo 

0,1 

bo 
bi 



P = 0.42 



0.43 /3 = 0.44 /3 = 0.45 



0.77(1) 0.82(1) 0.827(5) 0.70(2) 

1.2(1) 1.27(2) 1.42(2) 1.66(2) 

0.54(2) 0.55(2) 0.55(2) 0.58(1) 

0.25(3) 0.26(2) 0.27(1) 0.281(6) 



TABLE I. Fitting the density to the power law ansatz ||4T 
and to (|42l. 



We point out that the data are quite well fitted by the 
ansatz 



plrr? = hap/m -\- 61 (/i/m)^ 



(42) 



The coefficients 60: &i are also listed in ta ble [H Figure [7] 
shows the curves from the fit according to ( |41[ ), called fit 
1, and according to (42 1 called fit 2. The fitting curves 



have been generated using the fit parameters obtained 
for (3 = 0.42. 

It is interesting to note that a free Fermi gas at 
small, but non-zero temperature would suggest the /i- 



depcndence shown in (42 1. In this case, the fitting coef- 
ficient 60 should be proportional to the temperature. In 
order to explore the potential correspondence between 
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FIG. 10. Phase diagram: Magnetic susceptibility x per total volume (left) and density (right). 



the 0(2) quantum model and a Fermi gas, further sim- 
ulations for different lattice sizes and aspect ratios are 
necessary. This is left to future studies. We also point 
out that in the case of a free gas of hardsphere bosons [32 
one expects a linear rise of the density with the chemical 
potential /ic- Note, however, that the coefficient bi 
should be negative in the latter case. 



Let us finally explore the magnetic susceptibility x be- 
low and above the onset value of the chemical potential. 
Figure[8]shows x normalised to its zero temperature value 
Xo in order to remove the renormalisation constant. The 
inlay of this figure shows that x/xo is indeed independent 
of the renormalisation point specified by /?. Although the 
density remains zero for /i < m, the chemical potential 
impacts on the properties of the theory as indicated by 
the increase x/xo with ^. For /i > m, the data for x/xo 
largely depend on the value /3. Note that the simulations 
are carried out for a fixed number of lattice points. A 
change of /3 therefore also implies a change of the lattice 
volume. This is indeed the cause for the apparent scaling 
violations: Figure [9] basically shows the magnetic suscep- 
tibility per volume, i.e., x/Xo/inT-^y)- The volume V is 
the physical volume and the mass gap is introduced to 
render this factor dimensionless. The factor xo is again 
necessary to remove the renormalisation constant. We 
observe good scaling for the magnetic susceptibility in 
units of the (physical) volume for a wide range of /3 val- 
ues. An exemption are the results for /3 = 0.45 which 
might be too close to the phase transition at /3c- We 
interpret these findings as follows: For ^ > m, matter 
starts populating the ground state. This alters the in- 
teraction between the spins such that the U(l) global 
symmetry of the model breaks spontaneously leading to 
a superfluid phase. 

We are finally in the position to calculate the phase 



diagram of the 2 + 1 dimensional quantum 0(2)-theory 
in the fi/m and T/m plane. We have carried simula- 
tions with (3 = 0.44 using a lattice size of 64 x Nt, 
where Nf = 2 . . . 32. We performed simulations for 
/ia — 0.01 . . . 0.6 in steps of 0.01. Naturally, the temper- 
ature steps from (37) are not equally spaced. We used a 



cubic spline interpolation to generate a regularly spaced 
60 X 60 data set before plotting. 

To detect the superfluid phase in the phase diagram, 
we plotted the renormalised magnetic susceptibility in 
units of the physical volume, i.e., x/Xo/("i'^V^) where V 
is the (physical) volume of the 3-torus, i.e., V = L'^/T. 
The result is shown in figure [TOj left panel. As expected, 
we find the superfluid phase in the cold, but dense region. 
Increasing temperature at a flxed chemical potential fi > 
TO, rapidly dissolves the superfluid phase. Also shown is 
the density in the /i-T-plane. At low temperatures, we 
observe the onset at /i > to. For /i « to, we observe an 
increase of the density with temperature, which is due to 
thermal excitations with more particles overcoming the 
mass gap than anti-particles. 



V. CONCLUSIONS 

Using existing simulations techniques such as the Wolff 
cluster algorithm and the worm algorithm, we have thor- 
oughly analysed the 3d 0{2) model in the continuum 
limit P ^ 13c. In this limit, the model represents the 
quantum 0(2) model in 2 -|- 1 dimensions. The mass gap 
TO takes over the role of the free (scale) parameter, while 
the lattice spacing a is fined tuned by the /3 dependent 
correlation length, ma{/3) = l/^(/3). 

We have analysed the zero density limit by calculat- 
ing the (renormalised) magnetic susceptibility as a tunc- 
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tion of the temperature. As usual, the temperature de- 
pends on the extent Nt of the torus in time direction, 
i.e., T/m = /Nt. We find good scaling over more 
than a factor of 4 in the lattice spacing. The magnetic 
susceptibility is rather independent of the temperature 
for T < m/2 while a rather steep descent is observed for 
T > m. 

Charge is defined by means of the Noether current 
of the global U{1) symmetry of the theory, and finite 
densities are introduced using the chemical potential /i. 
In a first step, we analysed the /j dependence of the 
(renormalised) magnetic susceptibility and the density at 
"zero" temperature (isotropic lattice). Simulations are 
carried out using the worm algorithm. The onset chem- 
ical potential /ic is found to agree with the mass gap m 
of the theory. At /i > /i^, the density smoothly rises with 
increasing /i ruling out Bose-Einstein condensation. In 
fact, the ^ dependence of the density is more in line with 
a free Fermi gas (at small temperatures) . More investiga- 



tions are needed to understand the theory close to onset 
in terms of an effective theory. This is left to future work. 

For /i > /ic, the magnetic susceptibility scales with the 
volume and signals the spontaneous breakdown of the 
U{1) symmetry and therefore superfiuidity. We finally 
present the results for the full fj,-T phase diagram. In the 
cold, but dense regime we find a superfiuid phase, which 
is rapidly resolved by increasing the temperature. 
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